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ABSTRACT 

A single very deep 10' field has been imaged to r ~ 26. There are 2682 galaxies with 
23 < r < 26 with a size signihcantly larger than the seeing disk in the held. After 
correcting for telescope aberrations, possible guiding errors, and signal degradation due 
to seeing, I hnd a polarization signal of amplitude |p| = 2.7%, with an uncertainty of 
1.2% (95% conhdence limit). For the 1773 galaxies with 23 < r < 25, the measured 
amplitude is 2.4% ± 1.2% (95% conhdence limit). 

If this observed polarization is due to gravitational lensing by large scale structure, 
cosmological inferences can be made. The polarization amplitude would indicate that 
Slo X (78 ~ 1. 

The amplitude probability distribution is a Rayleigh distribution with variance so 
the probability T’(l/2 < |p|/(7p < 2) = 75%, and T’(l/3 < |p|/(7p < 3) = 93.5%. The 
measurement of the polarization amplitude in a single held can be used to constrain 
cosmological models since (Tp oc SIq x (Ts. 

Key words: Cosmology: observations - gravitational lensing - large scale structure of 
the Universe 
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ll INTRODUCTION 

o 

A^central question of modern cosmology is the degree to 
galaxies trace the underlying mass distribution of the 
universe. Large-scale inhomogeneities in the mass distribu¬ 
tion have been inferred from the dynamics of galaxies (eg. 
Faber & Burstein 1988, Bertschinger & Dekel 1989). How¬ 
ever, the determination of peculiar motions of galaxies is 
fraught with difficulties such as sensitivity to small errors 
in distance indicators and inhomogeneous Malmquist bias. 
As a result, observations that seem to indicate large-scale 
coherent flows of galaxies remain controversial. 

Weak gravitational lensing of distant galaxies, equiva¬ 
lent to a systematic induced image ellipticity without mul¬ 
tiple imaging, is an alternative probe of the large-scale mass 
distribution of the universe (eg. Kristian 1967, Gunn 1967). 
In this case photons are used as a probe of the potential in¬ 
stead of galaxies. Essentially, when light rays from a distant 
galaxy pass by an overdense region in the mass distribution, 
the observed images are slightly elongated tangentially with 
respect to the center of the perturbation. Such an effect has 
been observed in a number of rich clusters of galaxies, eg. 
Tyson, Valdes & Wenk(f990), Kaiser & Squires (1993). Un¬ 
derdensities in the mass distribution produce images with 
radial elongation. Independent perturbations along the line 


of sight add stochastically (Blandford et al. 1991 [BSBV], 
Miralda-Escud(= 1991, Kaiser 1992, Villumsen 1995). 

Early searches for weak lensing outside the region of 
rich clusters resulted in unsurprising non-detections (Kris¬ 
tian 1967, Valdes, Tyson & Jarvis 1983). More recently, from 
a single deep CCD frame of a blank held Mould et al. (1994) 
[M-94] attempted to measure the mean image polarization, 
p, of 4363 galaxies with magnitudes 23 < r < 26. Inside a 
circular aperture of radius 4.8', p = O.Of ± O.Of was found 
by M-94. 

As noted in M-94, there exists in the data a non-circular 
point spread function (PSF), possibly due to a telescope 
guiding error. Such an effect will lead to an artihcially in¬ 
duced image ellipticity that can mimic a cosmological grav¬ 
itational lensing signature. The artihcial signal is, however, 
present in both the galaxies and the stars (which will not 
be gravitationally lensed) and will be of least importance 
for large galaxy images. It can, in principle, be removed by 
an extrapolation to inhnite image size and an attempt to 
do this was carried out in M-94. In this Letter a re-analysis 
of the M-94 data is presented, focusing on the removal of 
the PSF from the galaxy images and optimal weighting of 
the data to hnd the mean of the distribution of image po¬ 
larizations. Provided it is possible to estimate the PSF from 
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the unsaturated stellar images, a mean image polarization 
of |) = 2.7% ± 1.2% (95% confidence bounds) is found. 

The theory of the removal of the PSF to obtain the true 
image orientation is outlined in §2. In §3 the observational 
data and limits on classical aberrations induced by the tele¬ 
scope and imager are discussed. In §4 the results of §2 are 
applied to the M-94 data using an optimal weighting scheme 
based upon the observed distribution of image orientations. 
A discussion of the implications of the results of §4 is given 
in §5. 


2 IMAGE POLARISATION AND PSF 
REMOVAL 

In the notation of BSBV the complex orientation of an image 
(galaxy or star) is dehned as x = {I 22 — In — Uln)/{In + 
I 22 ), where Lj are the image second moments. The modulus 
of X is (d^ — }J)l{a? b^) where and b^ are the principal 

second moments of the intensity. To linear order the mod¬ 
ulus of X for an image is equal to its ellipticity as dehned 
conventionally. The phase of x is twice the position angle, 
<t>. 

The effect I seek is the mean polarization, a net el¬ 
lipticity in the galaxy images induced by weak lensing by 
large-scale structure. In the limit of mildly elliptical images, 
as will be assumed, the mean polarization is p = (x), where 
the brackets ( ) denote an average over an area of sky large 
enough to contain many galaxy images but smaller than 
the size of the large-scale structure being probed. Assuming 
p ~ 0.03, as would be typical of an unbiased CDM universe 
with Qo = 1 (eg. M-94), ~ 1000 galaxy images are needed 
to measure this level of image polarization in the absence of 
systematic errors (eg. BSBV). 

To evaluate (x), the distribution of the true image ori¬ 
entations must be known. In the limit of noise-free data, 
the observed x for an image, x°t i® equal to the true X) 
X*. For observational data the primary difficulties in deter¬ 
mining X* ^re effects of the atmosphere (eg. seeing) and 
trailing of the images due to imperfect guiding. For suffi¬ 
ciently large images the effects of guiding errors and seeing 
can be ignored. However, for typical images seeing dimin¬ 
ishes the eccentricity (“circularization”) and guiding errors 
induce artihcial eccentricity which can mimic a lensing sig¬ 
nal. Possible problems with matching of frames (f.e. linear 
and rotational offsets) can be characterized by an asymmet¬ 
ric PSF. Such effects can be modeled as a convolution of the 
true images with a PSF and here we outline the procedure 
for removing the effect of the PSF from images comparable 
to or larger than the isophotal areas of the unsaturated stel¬ 
lar images. Using the unweighted image second moments, I 
show that if we know the second moments of the PSF (as 
characterized by the unsaturated stellar images) the map¬ 
ping from X* to x° is straightforward. 

Consider the effect of the PSF on a single image. The 
light prohle of the true image is fi{xi,X 2 ) and the PSF is 
characterized by a function S{xi,X 2 )- No specihc assump¬ 
tions about the forms of p and S are made other than that 
they both have unit integral and zero dipole moment. From 
the convolution of p and S an expression for the true un¬ 
weighted image second moments, /F, in terms of the ob¬ 


served image second moments, 7° , and the second moments 
of the PSF, Ajj, can be derived. 

A photon that would have arrived at the detector at 
{xiyXn) in the absence of the PSF has a probability den¬ 
sity <S(Axi, AX 2 ) of arriving at {xn -\- Axi,X 2 -|- AX 2 ). So 
the photon that would have contributed XtXj to the second 
moments instead contributes the term 

XtXj f d^Ax{xt-\-Axt){xj-\-Axj)S{Axi,Ax2) 


= J (JAx {xiXj -\- XtAxj -\- XjAxt -\- AxtAxj) 

xS{Axi,Ax 2) = XtXj-\-Atj. ( 1 ) 

Integrating over the entire image the result is 

It] = f d'^x (x^Xj + A,j) ii{xi,X2) = ijj + A,j. (2) 


That is, the observed image second moments are simply a 
sum of the true second moments (f.e. the moments one would 
measure in the absence of the effect of the PSF) and the 
second moments of the PSF. From this the observed image 
orientation, x°) can be parametrized in terms of the true 
image orientation, x*) and the second moments of the PSF: 

o _ I 22 ~ Ill ~ liln 

N Jo I Jo 
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I dehne a variable T = 1 — (An -|- A22)/(7n + I 22 ) for 
an image and the complex orientation of the PSF as x^ = 
(A 22 — All — 2fAi2)/(A22 + All). With this parametriza- 
tion the observed image orientation, x°t i® a simple linear 
function of T, x* j and x^ 

X°(T) = x^ + (x‘ - T. (4) 

Thus for T Si 0 (f.e. a very small galaxy) the observed x for 
an image is just that of the PSF while for a very large galaxy, 
T Si 1, X is unaffected by the PSF and therefore x° = X*■ 
Furthermore, the observed x is a linear function of T. This 
result follows from simple parametrizations and does not 
require any special assumptions about the image and PSF 
prohles. Additionally, it is not necessary to assume the PSF 
results from seeing and guiding errors. Neither is it necessary 
to assume galaxies follow a standard prohle. In principle eqn. 
4 allows us to infer the true x for a single image; however 
due to the presence of noise this is an unstable procedure 
and will not be attempted. 

For given values of T, I average x over the appropriate 
images and hnd the mean image polarization as a function 
of T is 

piT) = x'^ + (p‘ - X^) T. (5) 

Thus, if I measure the image polarization as a function of T 
and I also know x^j then the true mean polarization of the 
galaxies, p*, is the linear extrapolation of p(T) to T = 1. 

In observational data the brightness is not a continuous 
function, but rather the data is in the form of photons per 
pixel. This will induce an error in the second moments, how¬ 
ever, as long as an image is well sampled, as in the present 
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data, this error is negligible. A further problem is that the 
unweighted second moments are unstable to the presence of 
noise. To counter that instability, the data is usually cut at 
some surface brightness level and only contiguous pixels with 
a surface brightness above the critical level are counted in 
the second moments. This will lead to an error in the map¬ 
ping from true x to observed x (E<1- 4). For a given data set 
it is necessary, through Monte Carlo simulations, to measure 
how accurate this mapping is in the presence of noise. For 
the current data set, with a magnitude limit of r < 26, this 
mapping is accurate within the noise in the data and does 
not induce spurious polarizations. 


3 OBSERVATIONAL DATA 

The imaging data used are of a single 9.6' x 9.6' blank held 
centered on a(1950) = 17''21'"07" 5(1950) = -|-49°52'2l", 
taken in Gunn r. The data were acquired with the COSMIC 
imaging spectrograph (Dressier et al. 1995) on the 5-m Hale 
telescope during periods of good seeing (0.7"-0.9"). The h- 
nal stacked image consists of a total of 24.0 ksec integration, 
has a 1(7 surface brightness limit of /ir = 28.8 mag arcsec”^, 
seeing of 0.87" FWHM, and a total area of 90.1 arcmin^. 
An object catalogue was created from this frame using FO- 
CAS (Valdes 1982) and contains ~ 6600 objects brighter 
than the 80% completeness limit of r = 26.2. The reduction 
of the data to a catalogue of detected objects is detailed in 
M-94, which also estimates the mean redshift of galaxies in 
the range 23 < r < 26 to be 2 fa 3/4. 

An obvious source of spurious ellipticity in the images 
is classical aberration due to the primary mirror and im¬ 
ager. From several offset images of the globular star cluster 
M3, M-94 placed limits on the amount of classical aberra¬ 
tion in the data. The only signihcant aberration found was 
distortion (f.e. the classic “barrel” or “pincushion”), corre¬ 
sponding to a radial image displacement of-2.6" in the frame 
corners and —0.9" at the centers of the frame edges. Individ¬ 
ual frames were geometrically remapped prior to stacking, 
otherwise an artihcial polarization of ~ (At^/10'' f percent 
would have been created in the combined frame, where the 
frame to frame displacement is Ad ~ 10"-30". As in M- 
94, due to the signihcant remapping required in the frame 
corners, analysis of the data is limited to objects whose cen¬ 
troids are within a circular aperture of radius 4.8', centered 
on the chip. Results of direct measurements of the classical 
aberrations from the star helds agree well with a ray trace 
analysis of the COSMIC optical design using the ZEMAX 
code and the spurious polarization induced by all classical 
aberrations is estimated to be p ^ 0.01. The observed two 
point correlation function of the shear Cpp{9), see M-94, Fig. 
2, can be used to constrain the aberration induced polariza¬ 
tion. This induced polarization would lead to an anticorre¬ 
lation in Cpp{9) on large scale which is not observed. 

It is also possible for atmospheric effects to produce 
spurious ellipticities. From the width of the r hlter, M-94 
estimated that atmospheric refraction contributes a uniform 
polarization of at most 0.006 in the individual frames and 
the net induced ellipticity in the combined frame is negligi¬ 
ble. 

Sensitivity variations across the held due to inadequate 


hathelding can also induce spurious polarizations. However, 
the variations would have to be much larger than the esti¬ 
mated hathelding uncertainty to signihcantly inhuence the 
results. 


4 OBSERVED MEAN IMAGE POLARIZATION 

Here the results of §2 are applied to the M-94 data and an 
estimate of p* is obtained. The stars in the frame provide 
a measure of the PSF since they are intrinsically point-like 
objects that are not subject to gravitational lensing. The 
second moments of the unsaturated stellar images are cal¬ 
culated in the same manner as those of the galaxies and pro¬ 
vide measures of = ( —0.008±0.005)-|-(0.017±0.005)i and 
A 22 + All = 5.9 ± 0.1 (pixel units). These numbers are the 
crucial ingredients for the PSF removal. In a least squares 
ht, including linear and quadratic terms, to the PSF as a 
function of position, no statistically signihcant variations of 
the PSF across the held were found. 

Throughout the analysis all galaxy images with 23 < 
r < 26, isophotal area larger than the mean isophotal area 
of the unsaturated stars, and whose centroids are within 4.8' 
of the center of the chip are used. There are 2682 galaxies 
that satisfy these criteria. 

A single image can be seen as a measurement of the 
polarization with errors given by the intrinsic eccentricity 
of the source plus the noise associated with the measure¬ 
ment and 1 wish to hud the optimum estimator for the 
polarization. To do this 1 dehne p(T) = Re( {x{T)) ) -|- 
ilm((x(T))) = Px{T) -\- ipy{T). The 2682 galaxy images 
were divided into 10 equally-spaced bins in T and {xx{T)) 
and {Xy{T)), together with were used to determine 

Px{T) and Py{T) independently as follows. 

In order to use a maximum likelihood (ML) estimator of 
Px,y{T) the distribution / of Xx,y for a given value of T must 
be known. In Fig. 1 the distribution of \xy \ for galaxy images 
with 0 < T < 1 is shown. From this hgure, for T ^ 0.8 the 
distribution is nearly exponential with a turnover near zero 
and is flat for larger T. Similar results are obtained for the 
distribution of \xx\. The mean eccentricity is observed to 
be an increasing function of T. This is due to two effects, 
the larger images are less affected by the circularizing effects 
of seeing, and for a given image size the eccentricity is an 
increasing function of T. 

In an ML estimation of the polarization, the optimum 
weightfunction Wml is WmlO) ^ ~[0f!dx)!(f(x)-xy\.T}ae 
ML estimator then suggests weighting of the data to hud 
{Xx,y) where the weight function is of the form WMhlXx.y) = 
[1 + (Xx.y!^ • The uncertainty, 8, in individual mea¬ 
surements of Xx,y is estimated to be ~ 0.07. Such a weight¬ 
ing would indicate that the information is in the rounder 
images. However, the ML estimator may not be the opti¬ 
mal estimator since, in order to apply it, the second deriva¬ 
tive of the distribution function / must be well-behaved. 
This is not true for a peaked distribution as in Fig. 1, in 
which case it is better to search for the mode of the distri¬ 
bution. For this reason 1 parametrize a weighting function 
Wn(xx,y) = -\- (Xx.yOfY^^ and determine p(T) using 

several choices of n. More negative values of n correspond 
to a stronger weighting of the data towards round objects 
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Figure 1. Distribution of \xy \ in 10 bins of T. The ordinate is 
logarithmic. 

and n = 0 corresponds to no weighting. The mean and vari¬ 
ance of the Xx, 7 j distribntions are determined according to 
the standard dehnition of the weighted mean. For large val- 
nes of T the distribntions of Xxyj are no longer exponential 
and no weighting of the data is done to determine {Xx,y) in 
the largest T bin. 

From eqn. 5, the real and imaginary components of p*, 
Px,y, are estimated independently by standard weighted lin¬ 
ear least sqnares hts to the observed components of p(T) and 
X'®. An important point is that, when determining {xx,y(T)), 
the appropriate valne of x in FF„ is x* not x°■ Thns, best 
estimates of px,y are determined in an iterative manner. For 
a given valne of n, 1 begin by assnming x* = X°j nse x° ns x 
in W„ and determine hrst estimates olpx,y as the extrapola¬ 
tions to p{T = 1) of the least sqnares hts to the components 
of p{T). The hrst estimate of px,y is then nsed to estimate 
X* from eqn. 5, which is then nsed as x in W„, and new 
estimates of px,y are obtained nsing least sqnares and ex¬ 
trapolating to p{T = 1). From the new estimates of Px,y, 
new estimates of x* are obtained, and the entire procednre 
is repeated nntil convergence is reached (typically less than 
10 iterations). The procednre has, from Monte Carlo sim- 
nlations, proven to be stable and, from the data itself, has 
proven to be insensitive to the starting assnmptions oipx,y. 
In Fig. 2, the hnal hts to px,y{T) for the case of ra = —2 
and 8 = 0.07 are shown. The PSF data points are located at 
T = 0 and were inclnded in each iteration of the linear hts 
to p{T). The points at T = 1 indicate the inferred mean po¬ 
larization and the error bars are the formal estimates from 
the least sqnares ht. The straight lines in Fig. 2 are clearly 
good hts bnt they are, in reality, too good since the proce¬ 
dnre nsed to determine p{T) introdnces correlations between 
the data points, resnlting in an nnderestimate of the errors 
in Px,y To obtain more realistic estimates of the errors in 
Px,y, the galaxy and star images were bootstrap resampled 
1000 times and the corresponding best-ht px,y determined. 
Resnlts for n = (0, —1, —2, —3) are shown as scatter plots in 


Figure 2. Linear fits to px,y (T) for n = —2, 5 = 0.07. The open 
squares indicate galaxy data points, the open circles are the PSF 
data point. The hlled circles indicate the best ht values for the 
true polarization. All error bars are formal 1-a error bars. 
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Figure 3. Scatter plots of best ht true polarizationfor 1000 boot¬ 
strap samples of the data, 23 < r < 26, for n = (0, —1, —2, —3), 
S = 0.07. Large dot indicates best ht value of polarization. 

Fig. 3 and 1 take the variance in px,y to be the variance as 
determined from the resamplings. 

The optimal valne of n to nse in the weighting of the 
data to determine {Xx,y) has not yet been addressed. A nse- 
fnl criterion for the optimal valne of n is that which mini¬ 
mizes the variance in x*- From Fig. 3, the variance in px,y 
is smallest for n = —3. In Table 1 best estimates of the trne 
polarization and its associated error are listed for FF„ with 
n = 0, 1/2, —1, —4. For n < —2 the polarization and the 
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error are nearly independent of n and the formal best valne 
of the mean polarization is {px,Py ) = ( —0.5±0.5, —2.7±0.6). 
This error estimate inclndes the error indnced by the nncer- 
tainty in what is the optimal weighting scheme. When this 
analysis is repeated for the 1773 galaxies with 23 < r < 25, 
very similar resnlts are obtained, see Table 1. This indicates 
that the signal is not generated exclnsively by galaxies near 
the completeness limit. From Table 1 it is seen that the er¬ 
ror bars for the polarization determinations with different 
weighting schemes overlap. 

For r > 25 there are very few images signihcantly larger 
than the seeing disk. This means that most of the statistical 
signihcance of the polarization measnrement comes from the 
brighter images, f.e. 23 < r < 25, and that it is not possible 
to derive the trne polarization from the faintest sonrces, f.e. 
25 < r < 26, only. This indicates that going to very faint 
magnitndes is only meaningfnl in exceptional seeing condi¬ 
tions. 

23 < r < 26 23 < r < 25 


n 

Px 

Py 

Error 

Px 

Py 

Error 

0 

-1-0.4 

-1.6 

1.2 

-0.2 

-1.4 

1.2 

-0.5 

-1-0.3 

-1.7 

0.9 

-0.2 

-1.7 

0.9 

-1.0 

-1-0.2 

-1.9 

0.7 

-0.2 

-2.0 

0.8 

-1.5 

-1-0.0 

-2.2 

0.6 

-0.3 

-2.3 

0.7 

-2.0 

-0.3 

-2.5 

0.6 

-0.4 

-2.4 

0.6 

-2.5 

-0.5 

-2.7 

0.5 

-0.4 

-2.4 

0.6 

-3.0 

-0.5 

-2.7 

0.5 

-0.5 

-2.3 

0.6 

-3.5 

-0.5 

-2.6 

0.5 

-0.5 

-2.2 

0.5 

-4.0 

-0.5 

-2.5 

0.5 

-0.5 

-2.1 

0.5 


Table 1: Table of measnred “trne” polarizations {px,Py) (in 
%) with error estimate for different valnes of the weighting 
parameter n. Colnmns 2-4 are for the magnitnde range 23 < 
r < 26, while colnmns 5-7 are for the magnitnde range 23 < 
r < 25. 

The weighting scheme is nsefnl for hnding the mean 
polarization across the held. However, this same weighting 
scheme wonld lead to an artihcially hat Cpp{9) if applied nni- 
formly across the held. There are not enongh galaxies in the 
held to allow the measnrement of the variations of the polar¬ 
ization across the held nsing an optimal weighting scheme. 
Cpp{9) then has to be determined from the nnweighted data 
and the error bars are too large to make cosmological infer¬ 
ences. 


5 DISCUSSION 

I have measnred a statistically signihcant polarization in a 
single held. The most controversial aspect of the data redac¬ 
tion is the PSF removal. Under the assnmptions that system¬ 
atic errors have been dealt with properly and that the po¬ 
larization is indnced by gravity on cosmological scales 1 can 
make theoretical inferences. Since the variance of the polar¬ 
ization scales as the hrst moment of the power spectrnm of 
density hnctnations with a low pass hlter, the main contribn- 
tion to a polarization signal will occnr near the scale where 
the slope of the power spectrnm is -1. For a CDM, or HDM 
power spectrnm (Bardeen et al. 1986) the mean polarization 


for a 10* held, where the sonrces are at a typical redshift 
z Si 3/4, will be generated on scales, Tcdm = 27r/A;cDM ~ 
20Mpc/(QoA), and Thdm = 27r/A;HDM ~ 30Mpc/(QoA). 
Here h is the Hnbble constant in nnits of 100km s~^Mpc~^. 

Villnmsen (1995), §3.4, shows that the most likely po¬ 
larization amplitnde ||)| is ap, the sqnare root of the polar¬ 
ization variance, and that the probability P{l/2 < |p|/(Tp < 
2) = 75%, and P(l/3 < |p|/(Tp < 3) = 93.5%. Thns the 
measnrement of the polarization in a single held carries cos¬ 
mological information. From Villnmsen (1995), Fig. 1, we 
see that the implication in an HDM or CDM nniverse is 
that the prodnct of Qo and as is of order nnity. The resnlts 
in this paper shonld be conhrmed or related by other inves¬ 
tigations of this and other helds dne to the difhcnlt natnre 
of the observations and data rednction. 
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